Towards a measurement of quantum backflow 
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Quantum backflow is a classically forbidden effect consisting in a negative flux for states with 
negligible negative momentum components. It has never been observed experimentally so far. We 
derive a general relation that connects backflow with a critical value of the particle density, paving 
the way for the detection of backflow by a density measurement. To this end, we propose an 
explicit scheme with Bose-Einstein condensates, at reach with current experimental technologies. 
Remarkably, the application of a positive momentum kick, via a Bragg pulse, to a condensate with 
a positive velocity may cause a current flow in the negative direction. 

PACS numbers: 03.75.-b,67.85.-d,03.65.Ta 



Quantum backflow is a fascinating quantum interfer- 
ence effect consisting in a negative current density for 
quantum wave packets without negative momentum com- 
ponents [ij. Despite its intriguing nature - obviously 
counterintuitive from a classical viewpoint - quantum 
backflow has not yet received as much attention as other 
quantum effects. Firstly discovered by Allcock in 1969 
[l|, it only started to attract some attention in the mid 
90's. Bracken and Melloy i2| provided a bound for the 
maximal fraction of probability that can undergo back- 
flow. Then, additional bounds and analytic examples 
where discussed by Muga et al. Recently, Berry 

analyzed the statistics of backflow for random wave- 
functions, and Yearsley et al. 0] studied some specific 
cases, clarifying the maximal backflow limit. However, 
so far no experiments have been performed, and a clear 
program to carry out one is also missing. Two main dif- 
ficulties are the measurement of the current density (the 
existing proposals for local and direct measurements are 
rather idealized schemes [81 ) and the preparation of the 
necessary states. 

In this Letter we derive a general relation that con- 
nects the current and the particle density, allowing for 
the detection of backflow by a density measurement, and 
propose a scheme for its observation with Bose-Einstein 
condensates in harmonic traps, that could be easily im- 
plemented with current experimental technologies. In 
particular, we show that preparing a Bose-Einstein con- 
densate with positive-momentum components, and then 
further transferring a positive momentum kick to part of 
the atoms, causes under certain conditions, remarkably, 
a current flow in the negative direction. 

Let us start by considering a one dimensional Bose- 
Einstein condensate with a narrow momentum distribu- 
tion centered around hki > 0, with negligible negative 
components. Then, we apply a Bragg pulse that trans- 
fers a momentum fi.g > to part of the atoms Q , popu- 
lating a state of momentum hk2 = hki + hq (see Fig. [1]) . 
By indicating with Ai and A2 the amplitudes of the two 



FIG. 1. (Color online) a) A condensate is created in the 
ground state of a harmonic trap with frequency ui^; at t = 
we apply a magnetic gradient that shifts the trap by a distance 
d; b) the condensate starts to perform dipole oscillations in 
the trap; c) when it reaches a desired momentum hki the trap 
is switched off, and d) then the condensate is let expand for a 
time t; e) flnally, a Bragg pulse is applied in order to transfer 
part of the atoms to a state of momentum hk2. 



momentum states, the total wave function is 

-^{x, t) = t{j{x, t) {Al + A2 exp [iqx + iLp\) , (1) 

where we can assume Ai,A2 G M+ without loss of gen- 
erality (with Af + A2 ^ 1), being an arbitrary phase. 
All these parameters (except the phase, that will be ir- 
relevant in our scheme), can be controlled and measured 
in the experiment. Then, by writing the wave function 
of the initial wave packet as 

ipix^t) ^ (l){x,t)eTip[ie{x,t)], (2) 

the expression for the total current density, J<i,{x^t) = 
(h/m)ln\ [5'*V'I'] can be easily put in the form 

^J*(x,0 = {ve)p^ + [p* + \4>\'{Al - AD] , (3) 

with p^,{x,t) = \(l){x,t)\^ {Al + Al + 2AiA2Cos{qx + ip)) 
being the total density. Therefore, a negative flux. 
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J^,{x,t) < 0, corresponds to the following inequality for 
the density 

1 



sign[?7(x,i)]p*(x,i) < 



-\<j>{x,t)\\Ai^Ai), (4) 



\vix,t)\' 

where we have defined rj{x,t) = 1 + 2W9{x,t)/q. Later 
on we will show that r]{x,t) < corresponds to a clas- 
sical regime, whereas for ri(x, t) > the backflow is a 
purely quantum effect, without any classical counterpart. 
Therefore, in the quantum regime, backflow takes place 
when the density is below the following critical threshold: 

This is a fundamental relation that allows to detect back- 
flow by a density measurement. It applies to any class of 
wavepackets of the form ([T}, including the superposition 
of two plane waves discussed in 0, 0| • 

In order to propose a specific experimental implemen- 
tation, we consider a condensate in a three dimensional 
harmonic trap, with axial frequency uJx- We assume a 
tight radial confinement, u;± 3> lOx, so that the wave 
function can be factorized in a radial and axial compo- 
nents (in the non interacting case this factorization is 
exact) [10[. In the following we will focus on the ID axial 
dynamics, taking place in the waveguide provided by the 
transverse confinement, that is assumed to be always on. 
We define — ^JhjmuJx, a\_ = ^ h/muj^. 

The scheme proceeds as highlighted in FiglT] The con- 
densate is initially prepared in the ground state Vo of 
the harmonic trap. Then, at i = the trap is suddenly 
shifted spatially by d (FigH^) and the condensate starts 
to perform dipole oscillations (FigdJj). The correspond- 
ing wave function is given by a scaling transformation 

El 

%{:{x, t) =ipo{x- C(t)) exp [ik{t)x - i'y{t)] , (6) 

with ^(t) = d[l — cos{ujxt)], k{t) = {m/K)ujxdsnY{ujxt), 
the value of the phase 7(t) being irrelevant. Next, dXt = 
ti, when the condensate has reached a desired momentum 
mvi = hki = hk{ti), the trap is switched off (FigHJ:). 
For convenience we redefine the spatial coordinate as x — 
^{ti) — > X, and we set 71 = 7(ti). Therefore, the wave 
function reads 



ip{x, ti) = ipoix) exp {ikix — iji) . 



(7) 



At this point we let the condensate expand freely for a 
time t (Fig[TJl). Hereinafter we will consider explicitly 
two cases, namely a noninter acting condensate and the 
Thomas- Fermi (TF) limit fllj, that can both be treated 
analytically. In fact, in both cases the expansion can be 
expressed by a scaling transformation 

1 , ( X — V\t 



ipix,t) = 




(8) 



X exp 



'2fi''b{t) 



ikiX 1 




where h{t) represents the scaling parameter (for con- 
venience we have also redefined the time coordinate, 
t — ti — > t) and j3{t) is an irrelevant global phase. This ex- 
pression can be easily obtained by generalizing the scaling 
in [13, [l3l to the case of an initial velocity field. 

In the non interacting case, the initial wave func- 
tion is a minimum uncertainty Gaussian, ^^0(2;) — 
(l/7r3Y^)exp \—x^/{2a^)\, and the scaling parameter 
evolves as h{t) 



tribution we have ipo{x) = [(/i 



ni. Hi. For a TF dis- 
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for 

I a: I < Rtf = v2/V™^ ^-nd vanishing elsewhere, with 
giD — g3D/{2na\) [15j, and the chemical potential /i 
fixed by the normalization condition Jdx\^p\'^ = N, the 



latter being the number of atoms in the condensate 11 1 
In this case b{t) satisfies b{t) — uj'^/b'^{t), whose asymp- 
totic solution, for t » l/ojj:, is b{t) ~ \/2tuJx [l6| . 

Finally, we apply a Bragg pulse as discussed previously 
(FigH^). We may safely assume the duration of the pulse 
to be very short with respect to the other timescales of 
the problem fl^l- Then, the resulting wave function is 
that in Eq. ([T]), with the corresponding critical density 
for backflow in Eq. ([5]). We have 



</.(x) 



1 f ^ ^ (hki/m)t 



while the expression for the phase gradient is 

h b{t) \ b{t) ] 



(9) 



(10) 



that, in the asymptotic limit t 3> l/w^;, yields the same 
result V9 = mx/ht for both the non interacting and TF 
wave packets. Eventually, backflow can be probed by 
taking a snapshot of the interference pattern just after 
the Bragg pulse, measuring precisely its minimum, and 
comparing it to the critical density. 

Before proceeding to the quantum backflow, let us 
discuss the occurrence of a classical backflow. To this 
end it is sufficient to consider the flux of a single wave 
packet (before the Bragg pulse), namely J^(x,t) — 
{h/m)\4)\'^V9 . In this case the flux is negative for x < 
vi{t - b/b) x_{t) (see Eq. Then, by indicat- 

ing with i?o the initial half width of the wave packet 
and with RL{t) = —b{t)Ro + vit its left border at time 
t, a negative flux occurs when Rl < x-, that is for 
-Ro > vi/b{t). In the asymptotic limit, the latter rela- 
tion reads i?o > fvi/uJx (/ = 1, 1/v^ for the non inter- 
acting and TF cases, respectively). On the other hand, 
the momentum width of the wave packet is Ap H/Rq 
[igj , so that the negative momentum components can be 
safely neglected only when mvi ^ H/Rq. From these two 
conditions, we get that there is a negative flux (even in 
the absence of initial negative momenta) when i?o 3> a^; . 
This can be easily satisfied in the TF regime. In fact, in 
that case the backflow has a classical counterpart due to 
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FIG. 2. (Color online) Plot of the function F(q, A2) defined 
in Eq. (I13|l . For a given value of a, the maximal backflow is 
obtained for the value A2 that minimizes F. 
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FIG. 3. (Color online) Plot of the flux J* (a;). Backflow cor- 
responds to J* < 0. Positions are measured with respect to 
the wave packet center. 



the force F = —dx(QiD\Pib(x,t)\-^) implied by the repul- 
sive interparticle interactions 12| • These interactions are 
responsible for the appearance of negative momenta and 
backflow. Then, sufficient conditions for avoiding these 
classical effects are ki l/a^ and Rq < fvi/uj^- 

Let us now turn to the quantum backflow. In order 
to discuss the optimal setup for having backflow, it is 
convenient to consider the following expression for the 
current density. 



for A2 that solves dF/dA2\a = 0, that is for 



J^i,{x,t) = — \q(Al + A1A2 cos {qx + ip)) (11) 
m 

Ve [Al + Al + 2A1A2 cos {qx + If))] , 



that follows directly from Eqs. ([3]) and ([5]). Let us focus 
on its behavior around the center of wave packet, namely 
at X « (hki/m)t = dujxt (the following analysis extends 
to the whole packet if duJxt is much larger than the con- 
densate width). In the asymptotic limit, the phase gra- 
dient at the center is V0\c ~ fci, and the flux in Eq. pT|) 
turns out to be proportional to that of the superposition 
of two plane waves of momenta ki and k2 = ki + q. This 
limit is particularly useful because for two plane waves 
the probability density is a sinusoidal function and the 
critical density becomes a constant, which in practice 
makes irrelevant the value of the arbitrary phase we 
cannot control. Then, the condition for having backflow 
at the wave packet center is 

kiAj + k2Al + (fci + k2)AiA2 cos {qx + tp) <0. (12) 

Since all the parameters ki and Ai are positive, the mini- 
mal condition for having a negative flux is kiA\ + k2A\ < 
{ki + k2)AiA2 (for cos(-) = —1), that can be written as 



F(a, A2) = l + aAl-{2 + a)A2^Jl - ^1 < 0, (13) 

where we have deflned a = q/ki. The behavior of the 
function _F(q;, A2) in the region where F < is depicted 
in Fig. [5] In particular, for a given value of the relative 
momentum kick a, the minimal value of F is obtained 



2aA 



,^l-Al + {2 + a){2Al-l)=Q. (14) 



In order to maximize the effect of backflow and its detec- 
tion, one has to satisfy a number of constraints. In prin- 
ciple. Fig. [2] shows that the larger the value of a = q/ki, 
the larger the effect of backflow is. However, q cannot be 
arbitrarily large as it fixes the wavelength A = 2'K/q of 
the density modulations, which must be above the cur- 
rent experimental spatial resolution Ur, A ^ CTr, for al- 
lowing a clean experimental detection of backflow (see 
later on). In addition, as discussed before, fci should 
be sufficiently large for considering negligible the nega- 
tive momentum components of the initial wave packet, 
ki ^ l/flx- Therefore, since the maximal momentum 
that the condensate may acquire after a shift d of the trap 
is hki — muJxd, the latter condition reads d ^ ax- By 
combining the two conditions above, we get the hierarchy 
1 <C d/ax {2iT/a){ax/(Jr)- Furthermore, we recall that 
in the interacting case we must have Rq < fvi/uJx = fd 
in order to avoid classical effects. Therefore, given the 
value of the current imaging resolution, the non inter- 
acting case (i?o ~ o-x) appears more favorable than the 
TF one (where typically Rtf ^ cLx)- Nevertheless, the 
latter condition can be substantially softened if the mea- 
surement is performed at the wave packet center, away 
from the left tail where classical effects take place. 

As a specific example, here we consider the case of an 
almost noninteracting ^Li condensate [2^ prepared in the 
ground state of a trap with frequency lOx = 2tt x 1 Hz 
(yielding Qx — 38 /Ltm). Then, we shift the trap by 0? = 
80 fj,m, so that after a time ti = Tr/{2uJx) = 250 ms the 
condensate has reached its maximal velocity hki/m = 
LJxd ~ 0.5 mm/s. At this point the axial trap is switched 
off, and the condensate is let expand for a time t 3> 
Qx/'jJxd until it enters the asymptotic plane wave regime 
(here we use i = 1 s). Finally we apply a Bragg pulse 
of momentum hq — ahki, with a = 3, that transfers 
24% of the population to the state of momentum hk2, 
according to Eq. ^ {A2 = 0.49, Ai ~ 0.87) 9]. The 
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FIG. 4. (Color online) Density (solid) and critical density 
(dashed) for the case discussed in the text. Backflow occurs 
in the regions where the density is below the critical value, 
see Fig. O Positions are measured with respect to the wave 
packet center. 



resulting flux is shown in Fig. [31 where the backflow is 
evident, and more pronounced at the wave packet center. 
The corresponding density is displayed in Fig. HI where 
it is compared with the critical value of Eq. ([5]). The 
values obtained around the center are p!|,"" ~ 8% and 

Let us now discuss more thoroughly the implication 
of a finite imaging resolution ar- First we note that 
experimentally it is difficult to obtain a precise mea- 
surement of the absolute density, because of uncer- 
tainties in the calibration of the imaging setup. In- 
stead, measurements in which the densities at two dif- 
ferent points are compared are free from calibration er- 
rors and therefore are more precise. Owing to this, 
it is useful to normalize the total density p^i{x,t) — 
\4>{x, i)p [A\ + A\ + 2A1A2 cos {qx + ip)) to its maximal 
value pj^"^ ~ |0maa;|2 (^^2 + + 2A1A2) . In addition, 
we have to take into account that, due to the finite reso- 
lution ar [m , the sinusoidal term cos (qx + if) is reduced 
by a factor ^ = exp[— after the imaging. Then, 
by indicating with Xmin{t) the position of the density 
minima, we have 



[Xmin 1 1) 



„max 



A\ 



2(^1 A2 ■ 



where "exp" refers to the experimental conditions. In- 
stead, the normalized critical density is (close the wave 
packet center, where V9 « ki) 



^max 



\^ix,t)\^ 



Al 



q + 2ki 10™'^- 12 (A1+A2)' 



(15) 



so that, assuming that \(l){x,t)\'^ varies to be on a scale 
much larger than (7^ to be unaffected by the finite imag- 
ing resolution, the condition for observing a density drop 
below the critical value reads 



Al 



Al 



2CA1A2 a A1-A2 



Al 



Al + 2CA1A2 



a + 2Ai+A2 



(16) 



In particular, in the example case we have discussed, 
backflow could be clearly detected with an imaging res- 
olution of about 3 /im, which is within reach of current 
experimental setups. 

In conclusion, we have presented a feasible experi- 
mental scheme that could lead to the first observation 
of quantum backflow, namely the presence of a nega- 
tive flux for states with negligible negative momentum 
components. By using current technologies for ultracold 
atoms, we have discussed how to imprint backflow on a 
Bose-Einstein condensate and how to detect it by a usual 
density measurement. Remarkably, the presence of back- 
flow is directly signaled by the density dropping below a 
critical threshold. 
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